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Abstract — This paper studies cooperative schemes for the 
inter-cell interference control in orthogonal-frequency-division- 
multiple-access (OFDMA) cellular systems. The downlink trans- 
mission in a simplified two-cell system is examined, where 
both cells simultaneously access the same frequency band using 
OFDMA. The joint power and subcarrier allocation over the two 
cells is investigated for maximizing their sum throughput with 
both centralized and decentralized implementations. Particularly, 
the decentralized allocation is achieved via a new cooperative in- 
terference control approach, whereby the two cells independently 
implement resource allocation to maximize individual throughput 
in an iterative manner, subject to a set of mutual interference 
power constraints. Simulation results show that the proposed 
decentralized resource allocation schemes achieve the system 
throughput close to that by the centralized scheme, and provide 
substantial throughput gains over existing schemes. 

I. Introduction 

In traditional cellular networks, the base stations (BSs) in 
different cells independently control the transmission with 
their associated users. The inter-cell interference is avoided 
or minimized by adopting different frequency reuse patterns, 
which only allow non-adjacent cells to reuse the same fre- 
quency band. The frequency reuse factor is assigned to specify 
the rate at which the same frequency band can be used in 
the network. Due to emerging high-rate wireless multimedia 
applications, traditional cellular systems have been pushed 
towards their throughput limits. As a result, it has been 
proposed to increase the frequency reuse factor such that each 
cell can be assigned with more frequency bands to increase the 
attainable throughput. In the special case where all cells can 
share the same frequency band for simultaneous transmission, 
this corresponds to the factor-one or universal frequency 
reuse. However, with more flexible frequency reuse, the inter- 
cell interference control becomes an essential problem in 
cellular systems, which has recently drawn significant research 
attentions (see, e.g., QJ— E)). 

For multicell systems with a universal frequency reuse, 
two promising approaches have been proposed to resolve the 
inter-cell interference problem (see, e.g., JT] and the refer- 
ences therein): interference coordination and network MIMO 
{multiple-input multiple- output). In the former approach, the 
performance of a multicell system is optimized via joint 
resource allocation among all cells, based on their shared 
channel state information (CSI) of all direct and interfering 
links across different cells. Furthermore, if the baseband signal 



synchronization among the BSs of different cells is available 
and the transmit messages of different cells are shared by 
their BSs, a more powerful cooperation can be achieved in 
the downlink via jointly encoding the transmit messages of all 
BSs. In this so-called network MIMO approach, the combined 
use of antennas at different BSs for joint signal transmission 
resembles the conventional single-cell multiantenna broadcast 
channel (BC) [1]. In this paper, the former interference co- 
ordination approach is adopted due to its relatively easier 
implementation in practical systems. 

More specifically, we study the inter-cell interference coor- 
dination for a two-cell OFDMA downlink system with univer- 
sal frequency reuse. All BSs and user terminals are assumed to 
be each equipped with a single antenna, and thus the system of 
interest can be modeled as a parallel interfering SISO (single- 
input single-output) BC. Promising applications of this two- 
cell system model are illustrated in Fig. Q] which shows a 
geographically symmetric setup with two adjacent macrocells, 
as well as a non-symmetric setup with one macrocell and 
one inside femtocell. This paper investigates the joint power 
and subcarrier allocation over the two cells to maximize 
their sum throughput, for both centralized and decentralized 
implementations. Specifically, for the centralized allocation, 
with the assumption of a global knowledge of all channels in 
the network, we propose a scheme to jointly optimize power 
and subcarrier allocation over the two cells by applying the 
Lagrange duality method from convex optimization [5|. This 
centralized scheme provides a performance benchmark for the 
decentralized schemes studied subsequently. 

For the decentralized resource allocation, this paper pro- 
poses a new cooperative interference control approach, 
whereby the two cells independently optimize resource al- 
location to maximize individual throughput subject to a set 
of preassigned mutual interference power constraints, in an 
iterative manner until the resource allocation in both cells 
converges. Two types of interference power constraints are 
further examined: one is to constrain the total interference 
power across all subcarriers from each cell to the active users 
in its adjacent cell, termed joint subcarrier protection (JSP); 
and the other is to limit the interference power over each 
individual subcarrier, termed individual subcarrier protection 
(ISP). Also, the optimal resource allocation rules for each cell 
to maximize individual throughput with JSP or ISP are derived. 

The rest of this paper is organized as follows. Section II 
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Fig. 1. System model for two-cell applications. 

introduces the two-cell downlink OFDMA system, and formu- 
lates the optimization problem for resource allocation. Section 

III presents the centralized resource allocation scheme. Section 

IV proposes two decentralized schemes via the cooperative 
interference control approach with JSP and ISP, respectively. 
Section V presents simulation results and pertinent discus- 
sions. Finally, Section VI concludes the paper. 

II. System Model And Problem Formulation 

As shown in Fig. Q] we consider a two-cell system sharing 
the same frequency band with each cell having a downlink 
OFDMA transmission. We use m £ {1,2} to denote each 
of the two cells, which are referred to as the 1st and 2nd 
cells in this paper, respectively. For convenience, let the 1st 
cell refer to the macrocell and the 2nd cell refer to either 
the macrocell or the femtocell in Fig. Q] The total system 
bandwidth shared by the two cells is assumed to be B Hz, 
which is equally divided into N subcarriers (SCs) indexed 
by n £ A= {1, 2, N}. Each SC is assumed to be used 
by at most one user inside each cell and could be shared 
between two users individually selected from the two cells. 
In addition, the users in the network are indexed by fci £ 
Ai={l,2,...,ifi}in the 1st cell and k 2 £ A 2 = {1, 2, K 2 } 
in the 2nd cell, where K\ and K 2 are the total numbers of 
users in each corresponding cell. 

Furthermore, we denote the channel power gains (amplitude 
squares) from the two BSs to their respective users, saying 
users ki,k 2 , in each cell as h nkl and h n k 2 , respectively. 
The inter-cell interference channel gain from BSi to k% is 
denoted by g n k 2 , while that from BS2 to ki is by g n k x - We 
assume that the noise at each user's receiver has independent 
circularly symmetric complex Gaussian (CSCG) distribution 
over SCs with zero mean and variance a 2 = zqB/N, denoted 
by CA/"(0,C7 2 ), where z is the noise power spectral density. 
In addition, the transmit power allocated to user fci at SC n is 
denoted by p nkl - Thus, over all users and SCs in the 1st cell, 
we can define a power allocation matrix Pi (Ki-by-N) with 
the non-negative elements denoted by p n kxi n £ A, fci £ Ai. 
Pi is assumed to satisfy an OFDMA-based power allocation 
(OPA), in which there exists at most one element in each 
column being larger than zero and all the other elements 
are equal to zero. This OPA constraint can be expressed 
as Pi £ S x = {Pi >:0| PnkxPnk', =0,VAi ^A'i,Vn}. 



Similarly, we can define the power allocation matrix for the 
2nd cell as P2 £ S 2 (i^-by-iV matrix) under a similar OPA 
constraint. The nth columns of Pi and P 2 are denoted by two 
vectors pi„ £ Si n and p 2 „ £ S 2n , where St n (S 2n ) is drawn 
from the nth column of Si (S 2 ). 

With the above system model and assuming that the inter- 
cell interference is treated as additional Gaussian noise at 
each user's receiver, the signal-to-interference-plus-noise-ratio 
(SINR) of user fci at SC n in the 1st cell is given by 



nr.rn Pnk 1 h n k l 

SllSR nkl - ■ — -. 



(1) 



Similarly, the SINR of user k 2 at SC n in the 2nd cell is 
denoted by SINR nk2 ■ Thus, the achievable sum-rate of user 
k m £ A m , m £ {1, 2} is given by 



r k m = 



(2) 



We consider the weighted-sum-rate (WSR) in each cell i.e., 



Rr, 



^k m = 



Wk m r kml m £ {1,2}, 



(3) 



where w km is the (non-negative) rate weight of user k m in 
the mth cell. With individual transmission power constraint at 
each BS, the following optimization problem can be formu- 
lated to maximize the system throughput defined as 



max Ri + R 2 
Pi,p 2 



(4) 



subject to 



N 



Pnk m <P,f,m£ {1,2}, (5) 



^k m = l Z-^ n =l ' 

where P® s is the given power constraint at BS m , and 

P m eS m ,me{l,2} (6) 
is the OPA constraint for the mth cell. 

III. Centralized Allocation 

In this section, we study the centralized optimization for 
jointly allocating resources in the two cells so as to maximize 
the system throughput, which corresponds to solving Problem 
(0]l globally with constraints (0 and ©. For the centralized 
allocation, it is assumed that all channel gains in the network 
are collected by a central controller, which is capable of 
performing a centralized resource allocation and informing the 
allocation results to each cell for data transmission. 

Due to the non-convex OPA constraint and the non-concave 
objective function over Pi and P 2 , the optimization problem 
in (0]i is non-convex and thus cannot be solved efficiently 
for the global optimum. Nevertheless, the Lagrange duality 
method can be applied to this problem to obtain a sub- 
optimal solution. Interestingly, according to J6), it has been 
shown that a so-called "time-sharing" condition usually holds 
for resource allocation problems in OFDMA, and the duality 
gap for such problems solved by the Lagrange duality method 
becomes asymptotically zero as the number of subcarriers in 



the system becomes large. Accordingly, in the sequel, we apply 
the Lagrange duality method to solve Problem d4)Q 

First, we express the partial Lagrangian of Problem © as 



In, P2n, Ai,A 2 ) n . 
+ A 1 P 1 BS + A 2 P 2 BS , (/) 



where, for each SC n G A, 



Ln (Pin, P2n, \\, A 2 ) = Z)fc=l r "fcl 



Z)fe 2 =l W k2 r nk2 ~ Ai X^LiPnfci ~ A 2 Efe 2 2 =iPn 



A'l 



^2 



(8) 



and Ai, A 2 are non-negative dual variables associated with the 
power constraints in (0 with m = 1 and 2, respectively. The 
Lagrange dual function is then given by 

g(\ 1 ,\ 2 )= max L(P 1 ,P 2) Ai,A 2 ). (9) 

F*i £Si ,P 2 ££2 



Hence, the dual problem can be defined as 

min #(Ai,A 2 ). 

Ai>0,A 2 >0 

For a given pairs of Ai and A 2 , we have 



(10) 



'(Ai,A 2 ) 



ff„(Ai,A 2 ) + A 1 P 1 BS + A 2 P 2 BB , (11) 



where g n (Ai, A 2 ) , n £ A, is obtained by solving the following 
per-SC maximization problem 

5„(Ai,A 2 )= max L„ (pi n , p 2n , Ai, A 2 ) . (12) 

The maximization problem in (O is thus decoupled into 
N per-SC resource allocation problems given by (TlZi . Due 
to the OPA constraints, for one particular SC n, it can be 
simultaneously assigned to one pair of users (fci,fc 2 ) from 
the two cells when the resultant L n (pi„, p 2n , Ai, A 2 ) in © 
attains its maximum value (with the optimized p nkl and p nk2 )- 
This user pair can be obtained by searching over all possible 
combinations from users k\ £ Ai, fc 2 G A 2 . Thus, the optimal 
SC and power allocation that solves the problem in (fT2l is 



(fci,fc 2 )= argmax <^ max L^ lM) \, (13) 
fcieAi,fc 2 eA 2 Lp n fci>o,p„fc 2 >o J 

where (fci, fc 2 ) is the selected user pair to share SC n, and 

l(kiM) _ Wki r nkl + w k2 r nk2 - Xip nkl - A 2 p nfc2 (14) 

is obtained from (HJ. For a given pair of (k\, fc 2 ), the optimal 



Pnfcj and p„fc 2 to maximize L 



(ki,k 2 ) 



in (fl4] > have no closed- 



form solutions due to the non-convexity of this problem. 
However, an iterative search based on, e.g., Newton's method 
10 can be utilized to find a pair of local optimal solutions 
for p nkl and p nk2 . Then, we can check all possible user 
combinations to determine the optimal SC allocation according 
to ([P31 with optimized power allocation. 

After solving the per-SC problems in (TTZt for all n's, a 
subgradient-based method, e.g., the ellipsoid method, can be 

'Based on numerical results, the duality gap for our problem at hand is 
nonzero only for a negligibly small portion of the total number of randomly 
generated channels, even with a not-so-large number of subcarriers. 



adopted to solve the dual problem in dTob so that the power 
constraints in (0 at both BSs are satisfied. The details are 
thus omitted for brevity. Note that Problem (2) can be solved 
in polynomial time with an overall complexity with order 
O (Itr out N KiK^ltrin) in its dual domain. Specifically, for 
one particular SC, we search for K\Ki combinations of user 
pairs and determine the power allocation for each user pair 
with Itri n iterations. In addition, Itr out is the number of 
iterations for solving the dual problem in (fTOt . However, this 
centralized allocation needs a system level coordination with 
all channel conditions in the two cells, which is a demanding 
requirement for practical applications. In the next section, we 
propose decentralized schemes for resource allocation, which 
can be implemented by each cell independently. 

IV. Decentralized Allocation 

In this section, a new cooperative interference control ap- 
proach is applied to design decentralized resource allocation 
schemes for the two-cell OFDMA downlink system. In this 
approach, each cell independently optimizes its resource al- 
location to maximize individual WSR under its BS's own 
transmit power constraint, as well as a set of newly imposed 
constraints to regulate the leakage interference power levels to 
the active users in its adjacent cell. The above operation iter- 
ates between the two cells, until both cells obtain a converged 
resource allocation under their mutual inter-cell interferences. 
Specifically, two decentralized allocation schemes are studied 
in this section corresponding to two different types of inter- 
ference power constraints, namely JSP and ISP. 

A. Joint Subcarrier Protection (JSP) 

In this subsection, we solve the optimal resource allocation 
problem of maximizing the 1st cell's WSR subject to its BS's 
power constraint and a given JSP constraint to the active users 
in the 2nd cell. Similar problem formulation and solution apply 
to the resource allocation in the 2nd cell and are thus omitted. 

Consider the resource allocation problem in the 1st cell 
subject to the leakage interference constraint for the 2nd cell. 
In order to characterize the leakage interference to the 2nd 
cell, BSi needs to know the interference channel gains from 
it to all active users over different SCs in the 2nd cell. Let 
fc 2 = 7r 2 £ A 2 denote the active user at SC n in the 2nd 
cell, with the corresponding interference channel gain from 
BSi to fc 2 being g nk ^. It is then assumed that g nk2 has been 
perfectly estimated by user fc 2 in the 2nd cell and fed back 
to BS 2 . After collecting g nk . 2 for all n's from its active users, 
BS 2 sends these channel gain values to BSi (via a backhaul 
link connecting these two BSs). Note that if a particular SC 
n is not used by any user in the 2nd cell, the corresponding 
interference channel gain g nk2 sent from BS 2 to BSi is set to 
be zero regardless of its actual value, so that this SC can be 
used by the 1st cell without any interference constraint. 

To maximize the WSR of the 1st cell, the following problem 
is formulated as 

max R\ (15) 

PieSi 



subject to 

X* i*** 1 ?*' (16) 

— 'fci— 1 * — ^n— 1 

where T 2 is the given JSP power constraint for protecting 
all the active users in the 2nd cell. Note that T 2 limits 
the interference power averaged over all the SCs; thus, the 
corresponding resource allocation scheme is refereed to as the 
Average scheme for convenience. 

We assume the non-negative dual variables associated with 
(fT6l l and (fTTT i are A, /1. Similarly as in the case of centralized 
allocation, for a given pair of A, /1, Problem ([T5T l can be 
decoupled into N per-SC problems in its dual domain, and 
the optimal allocation for SC n in the 1st cell is derived as 

/ci = arg max < max Ljt > , (18) 

feiGAi lf>r. fcl >0 1 ) 

where Z/£ = ^log 2 (l + ^Pn kl ) - (a + ^) p nfcl) 
with /Ji = p n h 2 9nki + 0-2 being the interference-plus-noise 
power at SC n. Equ. (fT8l means SC n should be assigned 
to the user, denoted by ki, giving the highest value of L n ^ 
with the optimized p n fe,. By letting dL x ^ /dpnk, be zero 
and considering non-negative power allocation, the power 
allocation in (fT8l should be optimized according to 

_ ( w kl 1 i& \ + 

Pnk ^{N\n2 X + ^ h nkl ) ■ (iy) 

Thus, (fT8l and (fT9l ) together provide the optimal resource 
allocation rules at all SCs with fixed A and /1. 

Then, the ellipsoid method can be adopted to iteratively 
search over A > and fi > so that the constraints in ( fTol l and 
( fTTI i are simultaneously satisfied. This algorithm also bears a 
linear complexity order i.e., O(ItrKiN), where Itr denotes 
the number of iterations for updating A and /1. Similarly as 
d!5t . the resource allocation problem for the 2nd cell can 
be formulated to maximize i? 2 subject to the transmit power 
constraint P 2 of BS2 and the JSP constraint T\ to protect 
the active users in the 1st cell. For a given pair of T\ > and 
T 2 > 00 the per-cell resource allocation described above can 
be iteratively implemented between BSi and BS2 until the SC 
and power allocation in both cells converges. 

B. Individual Subcarrier Protection (ISP) 

In this subsection, we study the decentralized resource 
allocation with ISP. Similarly to the previous case of JSP, 
we merely present the solution to the optimization problem 
for the 1st cel. With the same objective function as (fT5T l 
and BS transmit power constraint as ([Tol . we formulate the 
current problem via replacing the JSP constraint in (IPTi by the 
following ISP constraint over each individual SC: 

Eti tai ^- T2>eA ' (20) 

2 The methods for properly setting Ti and T2 can be found in the journal 
version of this paper Q. 



where T 2 n is the interference power constraint for protecting 
the active user at SC n in the 2nd cell. 

Again, we apply the Lagrange duality to solve the per- 
cell resource allocation problem with ISP. Following a similar 
procedure as in JSP, we can derive the following optimal SC 
and power allocation rules: 

fci=argmax< max L n] ., >, (21) 

k%A, \0< PnH <T«/g n%2 nki y 

[( w kl 1 I21 \ + T£ 1 

p nkl = mm i[ - — — , I , (22) 

\ In 2 A h nk J g nk2 J 

where L^ ki = ^-log 2 ^1 + ^^Pnk^j -^Pnki, with A being 
the non-negative dual variable associated with d20l >. 

According to (l2TT i and (l22l . the optimal SC and power 
allocation can be determined for all n's with any given A > 0. 
Then, the bisection method [5] can be used to adjust A so that 
the BS transmit power constraint ([Tol is satisfied. 

Nevertheless, it is not computationally efficient to individu- 
ally optimize T£ (T™) for each SC, thus two special schemes 
are further identified. One scheme is to set T" = TJ l = 
+00, Vn 6 A, which means that each cell is not aware of its 
interference to the adjacent cell, named as the No Protection 
scheme. The other scheme is to set uniform peak interference 
power constraints over all SCs, i.e., T™ = T£ , T 2 " = T 2 * , Vn 6 
A, named as the Peak scheme. 

V. Simulation Results 

In this section, simulation results are presented to evaluate 
the performance of the proposed schemes for the two-cell 
downlink OFDMA system. It is assumed that B = 100 MHz 
and N = 64. In addition, all users' rate weights are assumed 
to be one, and the noise power spectral density zq is set 
to be —100 dBm/Hz. Assuming independent (time-domain) 
Rayleigh fading with six independent, equal-energy multipath 
taps, the frequency-domain channel gains {h nkl }, {h n k 2 }, 
{Snfci}' and {g n k 2 } are modeled as independent CSCG ran- 
dom variables distributed as CA/"(0,a), CAf(0,b), CAf(0,c) 
and CAf(0,d), respectively. For convenience, we normalize 
o = l, and adjust b, c and d to generate different channel 
models. Figs. |2]and[3] show the results for two macrocells with 
K\ = K2 = 8 and Fig. [4] for the case with one macrocell and 
one femtocell with K\ = 8 and A' 2 = 2 (cf. Fig. 1). 

Fig. [2] shows the system throughput, R\ + i? 2 , achieved 
by different interference power constraints T\ and T 2 using 
the proposed decentralized scheme with JSP (i.e., the Average 
scheme) in Section IV.A for one particular channel realization. 
The channel gains are obtained by setting 6=1 and c = d = 
0.2, while the transmit power limits at two BSs are set equally 
to be P BS = 1 watt. In this figure, we have marked one local 
maximum point obtained by the iterative search method in 
IfTl . Also, we have marked the system throughput obtained 
by the centralized scheme proposed in Section III (Optimal 
scheme). It is observed that the system throughput achieved by 
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Fig. 2. System throughput of the decentralized scheme with JSP: One random 
channel realization. 

the decentralized Average scheme is suboptimal as compared 
to that by the centralized Optimal scheme. 

Fig. [3] shows the system throughput against the average 
inter-cell interference channel gain for various schemes. The 
channels are generated via b — 1 and c = d = g, with g being 
the average interference channel gain ranging from 10~ 4 to 
1. The proposed decentralized Average scheme achieves the 
system throughput close to that by the centralized Optimal 
scheme for all values of g, when the searched optimized 
values of T\ and T2 are applied Q. If instead the preassigned 
values for T\ and T2 are applied, throughput degradations are 
observed to be negligible in the case of T\ = T2 = O.lp for 
the low inter-cell interference regime with small values of g, 
and in the case of T\ = T2 = O.Olp for the high inter-cell 
interference regime with large values of g. In addition, the 
Half scheme (each cell orthogonally uses half of the overall 
frequency band) and the No Protection scheme are observed 
to perform poorly for small and large values of g. Moreover, 
the Average scheme with JSP performs superior over the Peak 
scheme with ISP, especially when g becomes large. 

Finally, Fig. [4] shows the system throughput for a macrocell 
with a femtocell inside it. The channel gains are 6 = 5, 
c = 0.1, and d = 0.5. The transmit power constraint at 
the macrocell's BS is assumed to be 1 watt, while that 
at the femtocell's BS is changed from 0.02 to 2 watts. It 
is observed that all proposed centralized and decentralized 
resource allocation schemes outperform the No Protection 
scheme in the achievable system throughput, which eventually 
becomes saturated with the increased inter-cell interference. At 
low femtocell SNR, there exists a noticeable throughput gap 
between the Average and Peak schemes, which is due to the 
fact that when the femtocell suffers detrimental interference 
from the macrocell, the Average scheme can opportunistically 
allocate the femtocell transmit power to a small portion of 
SCs with best channel conditions. On the other hand, at 
high femtocell SNR, both Average and Peak schemes tend 
to perform close to the Optimal scheme. 

VI. Conclusion 

In this paper, the downlink cooperative interference control 
in a two-cell OFDMA system is investigated with central- 



Fig. 3. System throughput vs. inter-cell interference channel gain. 
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Fig. 4. System throughput vs. average femtocell SNR. 

ized and decentralized implementations for joint power and 
subcarrier allocation to maximize the system throughput. It 
is shown that the proposed decentralized recourse allocation 
schemes via the new approach of inter-cell interference power 
protection achieve a performance close to that of the central- 
ized scheme in various system settings. In addition, the joint 
subcarrier protection (JSP) with average interference power 
constraint is shown to achieve a larger system throughput 
than the more stringent individual subcarrier protection (ISP) 
counterpart with peak interference power constraint. 

References 

[1] D. Gesbert, S. Hanly, H. Huang, S. Shamai Shitz, O. Simeone, and W. 
Yu, "Multi-cell MIMO cooperative networks: A new look at interference," 
IEEE J. Select. Areas Communications, vol. 28, no. 9, pp. 1380-1408, 
Dec. 2010. 

[2] L. Venturino, N. Prasad, and X. Wang, "Coordinated scheduling and 
power allocation in downlink multicell OFDMA networks," IEEE Trans. 
Vehicular Technology, vol. 58, no. 6, pp. 2835-2848, July 2009. 

[3] W. Yu, T. Kwon, and C. Shin, "Joint scheduling and dynamic power 
spectrum optimization for wireless multicell networks," in Proc. 44th 
Conf. Info. Science Sys. (CISS), Princeton, NJ, March 2010, pp. 1-6. 

[4] R. Zhang and S. Cui, "Cooperative interference management with MISO 
beamforming," IEEE Trans. Signal Processing, vol. 58, no. 10, pp. 5450- 
5458, Oct. 2010. 

[5] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge, U.K.: 

Cambridge Univ. Press, 2004. 
[6] W. Yu and R. Liu, "Dual methods for nonconvex spectrum optimization 

of multicarrier systems," IEEE Trans. Communications, vol. 54, no. 7, 

pp. 1310-1322, July 2006. 
[7] B. Da and R. Zhang, "Cooperative interference control for spectrum 

sharing in cellular OFDMA systems," submitted for publication. 



